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Abstract 



Highly dynamic magnetospheric perturbations in neutron star environments can naturally 
account for the features observed in Gamma-ray Burst spectra. The source distribution, how- 
ever, appears to be extraga lactic. Although noncatastrophic isotropic emission mechanisms 
may be ruled out on energetic and timing arguments, MHD processes can produce strongly 
anisotropic 7-rays with an observable flux out to distances of ~ 1 — 2 Gpc. Here we show 
that sheared Alfven waves propagating along open magnetospheric field lines at the poles of 
magnetized neutron stars transfer their energy dissipationally to the current sustaining the field 
misalignment and thereby focus their power into a spatial region ~ 1000 times smaller than 
that of the crustal disturbance. This produces a strong (observable) flux enhancement along 
certain directions. We apply this model to a source population of "turned-off" pulsars that 
have nonetheless retained their strong magnetic fields and have achieved alignment at a period 
of 5 seconds. 

Subject headings: acceleration of particles - cosmology: miscellaneous - galaxies: evolution 
- gamma rays: bursts - MHD - pulsars: general 
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1. Introduction 

The lack of a precise determination of a distance scale to Gamma-ray burst (GRB) 
sources has greatly hindered our theoretical understanding of these objects. Much of what 
we know about these bursts is based on inferences drawn from clues provided by their 
spectra, including: (1) Most events exhibit rapid variability, apparently on time scales 
shorter than the best available instrument temporal resolution, with one burst exhibiting 
structure on a 200 /is timescale (Bhat et al. 1992). This variability indicates that the 
sources must be very compact, characterized by a length scale no larger than ~ 50 — 100 
km. (2) Bursts can last anywhere from tens of milliseconds to as long as 900 seconds, 
and usually have a complex temporal structure. This complicated time dependence would 
seem to favor mechanisms that invoke highly dynamic perturbations in otherwise stable 
environments. (3) Typical GRBs emit a substantial fraction of their power at energies 
in excess of ~ 1 MeV, which suggests that nonthermal processes are responsible for the 
emission of the 7-rays. In this regard, the inverse- Compton scattering of soft photons by 
relativistic particles has been very successful in reproducing the observed spectra (e.g., 
Pozdnyakov, Sobol' & Sunyaev 1977; Canfield, Howard & Liang 1987; Melia & Fatuzzo 
1989; Ho & Epstein 1989; Melia 1990a,b). 

It is possible to account for these observations by invoking a model in which the bursts 
originate within the magnetosphere of strongly magnetized neutron stars. In addition 
to being very compact, these environments are subject to magnetic fluctuations on sub- 
millisecond timescales, and the constituent particles can be energized nonthermally via the 
induced electrostatic forces. 

As is well known, however, the neutron-star paradigm must be reconciled with the 
uniform, yet spatially truncated GRB distribution observed by the BATSE experiment on 
CGRO (Meegan et al. 1992). These observations seem to rule out nearby (i.e. Galac- 
tic) single population models, and have therefore led to renewed speculation that GRBs 
originate at cosmological redshifts. But a naive estimate of the burst energy required for 
such distant sources yields a value that is significantly larger than that which a neutron 
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star could reasonably produce, unless the event was catastrophic (e.g., the coalescence of 
a neutron-star binary, Narayan et al. 1992), which does not seem to be borne out by the 
time history of typical bursts. 

A resolution to this apparent conflict was recently proposed by Melia & Fatuzzo (1992, 
hereafter MF), in which sheared Alfven waves generated near the polar cap of strongly- 
magnetized neutron stars produce streams of relativistic particles that are focused by 
the underlying magnetospheric structure. These energetic charges upscatter the radio- 
frequency photons (emitted at larger radii) into 7"^ cones aligned with the underlying 
magnetic field lines, resulting in an enhanced 7-ray flux along preferred lines of sight. This 
anisotropic emission is such that a pulsar glitch releasing ~ 10^^ ergs of energy could be 
viewed as a GRB out to a distance of ^ 1 Gpc. A key assumption of this scenario is that 
the Alfvenic power can indeed be emitted anisotropically. We show in this Letter that 
the required focusing is a natural consequence of the dissipational properties of sheared 
Alfven waves whose shear lengthscales (s ^ 10 cm) are much smaller than the size of the 
region (^ 10^ cm) encompassing the overall Alfven wave fluctuation. As such, this work 
strengthens the case for a non-catastrophic, cosmological origin of GRBs, and supports 
our view that an improved understanding of the micro-physical processes in neutron-star 
environments can indirectly, though significantly, infiuence our study of galactic evolution 
out to a redshift in excess of 1 — 2 (Tamblyn & Melia 1993). 

2. Sheared Alfvenic Wave Dissipation 

The general theory of sheared Alfven waves (SAWs) has been developed elsewhere 
(Melia & Fatuzzo 1992; Fatuzzo & Melia 1993). Here we consider the global properties 
of their dissipation. We idealize the unperturbed polar cap region as a fully ionized, 
homogeneous plasma threaded by a uniform magnetic field Bq = Bqz. The sheared Alfven 
waves may therefore be described by magnetic perturbations of the form 

= Ba{y) exp{ikz — iu>t)x , (1) 

where Ba{y) is an odd function that characterizes the shear geometry. For our purposes 
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here, we take the shear profile to be 

(Bao S>y>s 

Ba{y) = l Baog{y) \y\<s , (2) 
[ -Bao -S<y<-s 
where g{y) is a continuous function that satisfies the condition ^(±s) = ±1, and where 
s <^ S {s being the shear lengthscale and S the lengthscale of the encompassing plane wave 
regions) . For convenience, we define r] as the ratio of the nonsheared to sheared surface 
areas {rj = S/s in the present geometry). 

It is clear from the form of Ba{y) and Ampere's Law that an electric field Eaz must 
exist inside the sheared region {\y\ < s) until a sufiiciently strong current Jg is produced 
parallel to the underlying magnetic field Bq. For these waves, the equilibrium Goldreich- 
Julian particle density uq is insufficient to support a current large enough to short out E^^. 
Charges must therefore be copiously stripped ofi^ the stellar surface, thereby inducing a 
charged particle fiow to give the required Jg. Since the Alfven speed is Va ^ c, the waves 
travel with a phase velocity = uj/k ^ c. In order for Jg and the encompassing magnetic 
shear to remain in phase, the particle fiow must be relativistic, and thus, have an average 
density 

ng^^^l.7x 10^ ( ( . (3) 

We note that the stripped particles escape from the system by fiowing out along the open 
magnetospheric field lines. 

In standard pulsar theory, radio emission results from the coherent motion of "bunches" 
of electrons streaming along open field lines with Lorentz factors 7 ~ 10^~^. As such, strong 
transient radio emission is expected to be a natural byproduct of sheared Alfven waves 
if similar particle energies are reached, and if this emission is produced with front-back 
symmetry along the local field-line direction, a large fraction of the overall radio fiux will 
naturally be funneled back onto the polar cap. Taking into account the coherent nature of 
the processes responsible for pulsar emission, we parametrize the fiux impinging onto the 
stellar surface by Fj. = ^r]~'^{ng/nc)'^ {Lc/7rRpJ, where Lq and nc are the Crab pulsar 
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luminosity and magnetospheric number density, respectively, and where Rpc is the radius of 
the open field line polar cap. Assuming a stellar radius of i?* = lO^cm, Rpc can be related 
to the pulsar period P via Rpc = 1.4 x 10^ cm (P/ls)-^/^. With Lc = lO^^ 

ergs s and 

nc = 10^3 cm~3, this yields 

= 2 X 10»ergs c.- s- « (^) ( ' (^) (il) , (4) 

where r]~^ is the sheared flow "fllling factor" . The parameter ^ encompasses both geometric 
and emission uncertainties, and as such is poorly known. We note that if ^ becomes too 
small 0.1 for the range of parameters considered here), the wave dissipation lengthscale 
due to field line annhilation becomes much larger than R^, and the SAW mechanism 
becomes inefficient at producing 7-rays (see the discussion after equation [7]). However, 
since the Crab pulsar is itself very inefficient at converting spin-down energy into radio 
emission compared to typical pulsars, and since we have made the conservative assumption 
that the (coherent) radio fiiix scales as ry"^ (i.e., the square of the total immher of particles), 
it is reasonable to assume that ^ 0.1 (see also the discussion after equation [11]). 

The presence of Fr results in a radiative drag on the relativistic current-carrying 
charges. By analogy with MHD phenomena, the current driving electric field (E = Eazz) 
must be generated within the shear at the expense of the magnetic wave energy. However, 
the sheared waves are distinguished from pure MHD fiuctuations for two important rea- 
sons. First, the charges which generate Jg are constrained to always move along the same 
Bo field lines, so that SAWs cannot easily change their initial structure. Second, the simple 
concept of Ohm's law is not valid for the relativistic fiow inside the shear. Indeed, once 
the particles become relativistic, the current quickly decouples from the driving electric 
field, and since the radiative drag increases rapidly with 7 (the particle Lorentz factor), 
one might expect that a mildly relativistic flow will be favored by the system. 

Though Eaz depends on the microphysics of the shear (including all the annihilation 
processes, such as the tearing mode instability), its value may be estimated with a relatively 
simple argument under the assumption that the annihilation time scale within the shear 
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is the shortest of the relevant time scales. The strength of the electric field is limited by 
the rate at which the oppositely-directed magnetic fluctuations are driven together by the 
large magnetic pressure gradients associated with SAWs. Since the Alfvenic field lines are 
strongly coupled to Bq via fliix freezing with the charged medium, this transfer of Alfvenic 
power into the shear is dictated by the diffusion rate within the resistive plasma in the 
region \y\ > s, where the resistance is provided primarily by e~ /radio photon scatterings 
that occur with a frequency ~ Uph ctmc c, in terms of the photon number density Uph and 
the magnetic Compton cross section ctmc- With eo the characteristic radio photon energy, 
each event imparts a momentum ~ eo/c to the electron (moving nonrelativistically with a 
velocity Ve ^ eo/wec), which must therefore interact with ~ {meVec/eo){c/ve) photons in 
order to suffer a significant deviation to its path. Thus, the electron defiection frequency 
is 



~ 4.3 X 10" s-i e ("^v 



B. 



aO 



2 / P 



ruec- V103/ VlQi^G/ VlOcm 

where ctmc ~ err/ 4: when Bao ~ -Bo ~ 10^^ G (Dermer 1990). 

With a conductivity a ^ uq ^ IrxigV^ and a diffusion time scale = ^'kcjS'^ jc? (the well 
known MHD value which is valid as long as ^ 27r/a;), the diffusion velocity Sjrj^ 
for the magnetic field lines is given as 



Vd =min 



c,8.7 X 10^ cm s"^ 
2 / B. 



(loO 



P 

5s 



Bo 



max 



5,3.5 X 10^ cm 
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aO 
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1012 Gy VSs/ VlO^^Gy VlOcm 
Thus, since ~ c, the waves dissipate over a lengthscale 



S 



104 



cm 



10 cm/ 



S 



104 



cm 



(6) 



(7) 



As long as i?^ ^ -R* (i-e., ^ is sufficiently large), most of the wave energy is channeled 
into the shear before the waves break, and we may equate the Alfvenic luminosity (~ 
B^QcAy^l^'K) generated at the stellar surface with the power (~ EAzJs-^sRd) dissipated by 
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the current as it converts magnetospheric energy into upscattered radiation. Here, and 
Ag are the surface areas corresponding to the wave and shear regions, respectively. This 
yields an average electric field strength 



E^g ftimin 



5 X lO^^sV cm-i ( -^^) ( ( ( , 1.4 x 10^° sV cm-^ 

^ ^ (to^) (toi2g) (si) (toi2g) (to^) (to^) ] ■ 

As before, we assume a typical pulsar spectrum specified as a steep power law with 
(fiiix density) index // above a break at frequency eo/h ~ 500 MHz. With 7 ^ 1, a 
lab frame photon with energy e will be blue-shifted to ~ 2je in the electron rest frame, 
which is well below the resonant energy = (So/44.14 x IQ-'^^G) mgC^, and its angle of 
propagation relative to the particle direction (and hence Bq) is ~ As such, e^/e ~ 
6x10^(^0/10^2 G)(e//i500Mhz)-i < 7^, and (Jmc ~ 4(jT(7e/eB)^ (Melia & Fatuzzo 1989; 
Dermer 1990). Balancing the accelerating force cEaz by the radiative drag ^'^FraMc/c, 
one therefore obtains 



1012 g; yiO^^Gj \5sJ VlO^cmy V/i500MHz 

(9) 

where /(^) = min[l, (C/Co)~^^^] and the parameter ^0 is defined as the smallest value of ^ 
for which Vd = c (e.g., ^0 = 35 for the represented parameter space), and where we have 
used T] = S/s. This result is consistent with the assumptions discussed above (e.g., that 
the particle motion is relativistic and sufiiciently energetic to produce the required radio 
luminosity and that (Jmc have an dependence). Ultimately, the current Js must decay 
in concert with the Alfvenic magnetic field, even though the particles remain relativistic. 
Evidently, the initially fully charge separated regions must merge together and neutralize. 
This behavior is expected since the power transferred from the wave to the particles is 
reduced as the magnetic field decays. As such, an increasing number of charges undergoing 
collisions will not be energized quickly enough to remain in phase with the wave, and are 
therefore swept up by the lagging oppositely charged wave region. 



The scaling of the magnetic fields Bq and Bao in the above equations was chosen for 
convenience and not to suggest that Bao ~ Bq. We note, however, that the presence of 
significant reconnection in the nonlinear regime would have the desirable effect of enhancing 
the annihilation rate within the shear (see the paragraph preceeding Eq. 5). 

3. The Cosmological Radio Pulsar Model For Gamma-ray Bursts 

In applying the above discussion to the cosmological gamma-ray burst model, we must 
now generalize to a more realistic magnetospheric geometry in which the field lines are 
more or less radial close to the stellar surface (cf. Eq. 7). Although globally the results of 
§2 are expected to apply here, we note an important difference between the two geometries. 
The underlying magnetic field strength Bq now decreases away from the stellar surface as 
(i?*/r)^, whereas in the absence of wave damping, the Alfvenic field decreases only as 
[R^/r), for which the waves eventually lose their linearity and shock. However, since the 
wave energy diffuses into the shear on a length scale < i?* , the Alfvenic field also decreases 
roughly as (i?*/r)^. 

It is evident from §2 that SAWs focus their energy into the internal current fiow. Thus, 

as long as 7"^ ^ s/R^, the flux is enhanced by a factor t] in certain directions, for which 
a source at a distance D will have an observable 7-ray fliix 

F, = rj^c (^X ^1.3x 10-^ ergs cm'^ s-UA,)( ' ( . (10) 

The observability of these bursts at cosmological distances imposes strict (but not unre- 
alistic) conditions on the model parameters, such as the required burst power [Lhurst ~ 
10^^ — lO'^^ ergs s~^), whose magnitude depends on whether the SAWs are generated only 
near the polar cap (where the field lines are most strongly coupled to toroidal crustal 
activity), or are generated throughout the entire stellar surface. 

Since the particle flow remains optically thin to the radio photons impinging upon the 
star, the dependence of the cross-section (see above) implies that the incipient radio 
spectrum is upscattered to a 7-ray spectrum with (power) index /x -|- 2 -|- 1, and very 
importantly, that the spectral radio break at eo is translated to the corresponding 7-ray 
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break at 

- - " ^'^^ ( (life) (1^) ( 10^) «K) ■ 

(11) 

independent of eo (we have again used 77 = ^/s). Here, ^(^) = min[l, (^/^o)""*^^^]- This 
result compares favorably with the observed value of ebreak (which after redshift is taken 
into account is seen to fall within the range ~ 100 keV —3 MeV; Schaefer et al. 1992), and 
suggests that ^/^o may be as large as 100. A more detailed description of the resulting 
7-ray spectrum is given in MF (see, for example. Figure 2 therein). For completeness, 
we note that a cylindrical shear would correspond to 77 ~ (S/s)'^ and the parameter S in 
Equation (11) should be replaced by s. Such a strictly confined shear region would thus 
appear to be unlikely, though it cannot be ruled out without a more detailed calculation. 

The question of which mechanisms contribute to the generation of SAWs is best ad- 
dressed by considering the constraints imposed on the model by the different probabilities 
of observing extragalactic and galactic events. If the source distribution is cosmological, 
the BATSE data imply a detection rate T^-c ~ 6 x 10~^ bursts per year per galaxy, with an 
actual rate TZt = T^c/Pci where Pc is the probability of detecting the extragalactic burst. 
Assuming that the number of galactic sources is roughly equal to the average number of 
sources in all other galaxies, we should therefore expect to see a galactic rate TZg — RrPg, 
where Pg is the corresponding probability of detecting a galactic burst in progress. Ac- 
cording to Equation (10), a galactic burst would be quite distinguishable, exhibiting fluxes 
of order 1 ergs cm~^s~^, i.e., ~ 10^ times larger than their extragalactic counterparts. 
Since none of these have ever been detected, we infer that ^ 0.01 per year, for which 
Pc/Pg > 6 X 10-^ 

What this means in practice is best seen with recourse to a specific scenario. Let 
us assume that the underlying sources are turned-off pulsars (i.e., neutron stars older 
than ~ 10^ years), which have nonetheless retained their strong magnetic fields and have 
achieved alignment at a period P ~ 5s, corresponding to an open field line polar cap radius 
Rpc Ri 6 X lO^cm. The stellar rotation implies that a given 7-ray flux region sweeps in and 
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out of the observer's view, but as long as the duration of an individual sweep is longer than 
the instrument resolution time r^, the inferred (average) flux is correctly given by Eq. (10). 
If, however, the sweep time is shorter than r^, the observed (average) flux for the most 
distant bursts drops below the instrument sensitivity. This means that an extragalactic 
burst is observable predominantly within a cone centered about the rotation axis, with an 
opening angle corresponding to a radius Rob = {s/T^){P/^i) (two poles). Thus, since the 
probability that a sheared region occurs within the area enclosed by this "observable" cap 
is ~ Rob/S (where Rob < S), the probability of seeing a burst in progress at cosmological 
distances must be 

The actual event rate would therefore be Rt ~ 10^ per galaxy per year, which in turn 
implies a stellar repetition time scale of ~ 10^ — 10^ years if the population is comprised 
of ~ 10^ - 10^ objects. 

The probability of detecting such a burst in progress within the galaxy is higher since 
for these events Rob ~ Rpc (see above). As such, Pg is simply the ratio of solid angles 
corresponding to the polar cap region and the entire star (i.e., An). For a P = 5s rotator, 
we therefore infer that a; 2 x 10~^, so that Pc/Pg ~ 4 x 10~^ for this population, which 
is consistent with the constraint discussed above. To put this result in another way, we 
would anticipate that such a population of GRB sources should produce an observable 
galactic "super" burst roughly once every ^ 50 years. In addition, we note that GRBs 
originating from extragalactic sources are expected to be accompanied by ~ 0.01 — 1.0 
Jansky radio bursts (see MF). For galactic bursts, the observable radio flux will be roughly 
10 orders of magnitude larger and therefore (as is the case for the 7-ray signal) quite 
distinguishable. Assuming that the radio emission is produced at ~ 1^R~^ above the polar 
cap, the probability of seeing the radio burst will be roughly 10 times greater than of seeing 
the corresponding 7-ray burst, suggesting an observable galactic "super" radio burst rate 
of roughly once every ^ 5 years. This result does not conflict with current observations 
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since only a small fraction of the sky is monitered by radio telescopes at any given time. 
A possible link between GRB sources and soft 7-ray repeater events, which in this picture 
would be interpreted as bursts viewed outside of the open field line cone, has been discussed 
elsewhere (Meha & Fatuzzo 1993). 

We are grateful to the anonymous referee for suggesting several improvements to the 
manuscript. This research was supported by NSF grant PHY 88-57218, and the NASA 
High-Energy Astrophysics Theory and Data Analysis Program at UA, and by the Compton 
GRO Fellowship Program at UM. 
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